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Agenda 


• Place-holder for agenda with times, etc 



Earth-GRAM Summary 

Earth-GRAM is a FORTRAN software package that can run on a variety of 
platforms including PC's 

For any time and location in the Earth s atmosphere, Earth-GRAM provides 
values of atmospheric quantities such as temperature, pressure, density, 
winds, constituents, etc 

Dispersions (perturbations) of these parameters are also provided and have 
realistic correlations, means, and variances - useful for Monte Carlo 
analysis 

Earth-GRAM is driven by observations including a troposphenc database 
available from the National Climatic Data Center 
Although Earth-GRAM can be run in a ‘stand-alone* mode, many users 
incorporate it into their trajectory codes 

The source code is distributed free-of-charge to eligible recipients 
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Background 


Original GRAM was developed (1974) to study Space Shuttle entry 
guidance and thermal protection issues 

GRAM is usually referred to as Earth-GRAM, to distinguish it from 
models of other planetary bodies (Venus. Mars. Titan, Neptune) 
Several Earth-GRAM updates have been made The current version 
is Earth-GRAM-99, version 3 (a new version is currently undergoing 
testing) 

Earth-GRAM99 is based on global climatologies from balloon, 
aircraft, rocket, ground-based remote sensing, satellite remote 
sensing, and in-situ satellite observations 

Earth-GRAM99 covers all global locations, all heights (surface to - 
1000 km), and all months 

Earth-GRAM 99 is diagnostic, not prognostic(i e it descnbes the 
atmosphere, but it does not forecasts) 
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Applications 

• Shuttle re-entry studies 

- Guidance algorithm design 

- Thermal Protection System design 

- RCS thruster fuel use analysis 

• Other NASA Projects (X33, X37, X38, X43, etc.) 

• Columbia investigation 

• Military applications 

• Stardust & Genesis missions 

• Constellation 


Available Documentation 


• The NASA/MSFC Global Reference Atmospheric Model - 
1999 Version ( GRAM99 ) NASA/TM-1 999-209630 [On the 
GRAM CD] 

• The NASA/MSFC Global Reference Atmospheric Model - 
1995 Version ( GRAM-95 ) NASAfTM4715 fOn the GRAM 
CD] 

• Bibliography [In Backup Slides] 

• WebPages: 

- http //see msfc nasa gov/tte/modei_gram htm 

- http //see msfc nasa gov/ModelDB/ModelDB htm 
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Earth-GRAM99 Output 


Earth-GRAM99 Output (cont’d) 

For any position and time, Earth-GRAM provides 


GRAM also provides mean values (no 

mean and (optional) perturbed values of: 


perturbations) of: 

• Temperature, K 


• Water vapor pressure 

• Pressure, N/m 2 


• Water vapor density 

• Density, kg/m 3 


• Relative humidity 

• E-W wind, m/s 


• 0 3 , N 2 0, CO, CH 4 , C0 2 , N 2 , 0 2 , O, Ar, He, H, N 

• N-S wind, m/s 



• Vertical wind, m/s 


GRAM does not provide global distributions of 



cloud cover, precipitation, visibility, or lightning 

7* *00X7 A' 
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Earth-GRAM99 Data Sources 



In addition to 
global data sets shown 
at left, Earth-GRAM99 
has an option to use 
Range Reference 
Atmospheres 
(1983 RRAs) at and 
near several launch 
range site locations 

Depending on site, 
RRAs cover either 
0-30 km or 0-70 km 
altitudes 


Source of GUACA Data - Daily Synoptic Observations 



• Darker Blue: 
standard 
rawinsonde 
sites with 
balloons every 
day, 00Z and 
12Z 

• Lighter Blue 
supplemental 
sites with 
balloons at 
irregular days 
and times 

• Satellite data 
used mostly 
over oceans 
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GUACA - Daily Data Gridded and Averaged 


GUACA Grid (CONUS Area Only) 



Longitude. dog. 


Gridded to 2.5-by- 
2.6 Deg Lat-Lon Grid 
(Used to initialize 
ECMWF weather 
forecast model) 

Data for GUACA 
assembled into 
monthly averages 
and standard 
deviations 
Data at Surface and 
14 pressure levels 
(1000 mb to 10 mb; 
- 0-27 km) 
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Other GUACA Data Characteristics 


• Binary data files containing gridded monthly 
means and standard deviations (sigmas) 

• Monthly means and sigmas for Period-of-Record 
(POR) 1980 - 1991 

• Monthly means and sigmas for individual years 
1985- 1991 

• Optional GGUAS (ASCII) data available for POR 
1980 - 1995 (no individual years) 

• Further description in files guaca.txt and 
gguas.txt (on the GRAM CD) 

7/20/2007 12 



Middle Atmosphere Program (MAP) 


Marshall Engineering Thermosphere (MET) 

• Height range: 20 km to 120 km 

• Monthly mean values of atmospheric 
parameters 

- 10° latitude interval from -90° to +90° 

- 20° longitude interval from 180° W to 160° E 

- 5 km height increments 

• Standard deviations depend on month, height, 
and latitude 

• Data from rockets and remote sensing 

• The dataset resides in the atmosdat.txt file 


• Used for heights above 90 km 

• Uses the Jacchia model to calculate 
atmospheric temperature and pressure 
based on solar angle/activity 

• Computes density from ideal gas equation 

• Computes winds using a geostrophic 
assumption with a viscous modification 

rawx? 


7/20/2007 


Range Reference Atmosphere Data (1983) 

• Site specific monthly (and annual average) data 

• Period of Record: 1957 - 1979 

• Included in GRAM software package 

• Sites: 

- Ascension Island 

- Barking Sands, Hawaii 

- Cape Canaveral, Florida 

- Dugway Proving Ground, UT 

- Edwards Air Force Base, CA 

- Eglin AFB, FL 

- Kwajalem Missile Range 

- Point Mugu, CA 

- Taquac, Guam 

- Vandenberg AFB, CA 

- Wallops Island, VA 

- White Sands. NM 


7/20/2007 


15 




Range Reference Atmosphere Data (e.g. tablel.wsm) 


JANUARY 

STATION - 722696 WHITE SANDS. NM 

Z MEAN U SD O R(U,V) MEAN V SD V MEAN WS SD WS SKEW WS NOBS 

0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

KM M/S M/S 

000 0.00 0.00 0 

246 0.02 2.29 - 

000 2.30 3.41 0 

000 6.61 5.60 0 

000 10.13 7.54 0 

000 13.01 9.96 0 

000 15.16 11.80 0 

000 16.62 13.52 0 

000 18.09 14.79 0 

000 19.64 15.73 0 

000 21.61 15.80 0 

000 23.37 15.51 0 

000 25.25 15.28 0 

000 25.54 13.44 0 

000 24 67 11.27 0 

000 22.62 10.52 0 

M/S M/S M/S M/S 

0000 0.00 0.00 0.00 0.00 0.00 

0457 0.06 2.77 2.73 2.33 1.28 1 

3109 -0.60 4.94 5.61 1.17 0.97 1 

1951 -1.61 5.75 9.66 4.15 0.45 3 

2764 -2.22 7.01 13.43 5.75 0.35 3 

2958 -2.13 8.56 16.94 7.73 0.59 3 

3858 -1.73 9.88 19.48 9.48 0.6S 1 

4683 -1.24 11.44 21.70 10.95 0.77 1 

5126 -1.15 13.05 23.89 12.09 0.74 1 

5153 -0.92 14.47 25.89 13.13 0.60 1 

4808 -0.87 15.88 27.83 13.91 0.50 2 

4369 -1.02 15.50 29.07 11.49 0.39 2 

3802 -0.61 14.47 30.20 12.95 0.41 2 

3053 -0.74 13.06 29.42 11.75 0.45 2 

4086 -0.92 11.77 27.81 10.07 0.12 2 

4217 -0.97 10.46 25.35 9.47 0.09 2 
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Range Reference Atmosphere Data (e.g. table2.wsm) 


JANUARY 

STATION - 722696 WHITE SANDS. NM 


Z MEAN P 

KM MB 

0.000 1021.6001 

1.000 903.8899 

1.246 876.9900 

2.000 800.3301 

3.000 707.1201 

4.000 623.2300 

5.000 547.8601 

6.000 479.9700 

7.000 418.9099 

8.000 364.2800 

9.000 3 IS . 2500 

10.000 271.5200 

11.000 232.9700 

12.000 199.3400 

13.000 170.1300 

14.000 145.0000 

15. 000 123.5600 


SD P SKEWP 
MB DEG 

9.0775 -0.24 280 
5 . 3868 -0.11 277 
4.8866 -0.04 276 
4.1375 -0.04 277 
4.1898 -0.26 272 
4.6670 -0.45 267 
5.1809 -0.S1 261 
5.4018 -0.51 254 
5.5981 -0.51 246 
5.6206 -0.45 239 
5.4136 -0.15 211 
5.1199 -0.27 224 
4.6069 -0.14 219 
3.8999 -0.04 215 
1.1878 0.0S 213 
2.5366 0.02 212 
1.9738 -0.05 210 


NT SD T SKBWT 

K DEG K 

.46 10.64 0.27 

.81 7.47 0.15 

.95 6.77 0.10 

.70 5.21 -0.98 

.66 5.00 -0.85 

.32 4.71 -0.41 

.09 4.54 -0.31 

.16 4.38 -0.30 

.71 4.33 -0.29 

.23 4.21 -0.15 

.76 3.79 0.01 

.80 3.24 0.01 

.19 3.16 0.32 

.43 4.76 0.47 

97 4.99 -0.05 

.47 4.03 -0.35 

.13 3.58 -0.33 


MEAN D 

G/M3 

1268.0000 

1132.0000 

1102.0000 
1003.0000 

902.6001 

811.7000 

730.7000 

657.7000 

591.3999 

530.3999 

471.8999 
420 8000 
370.0000 
322.5000 
277.2000 

237.9000 

204 . 9000 


SD D . . 
G/M3 
58. 9000 
35.4300 
30.8800 
20.1900 
14.4900 
10.4100 
8.0810 
6.4520 
5.3440 
4.8300 
4.9990 
6.3920 

8 5270 
10.4400 

9 6650 
7.1000 
5.8340 
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Range Reference Atmosphere Data (e.g. table3.wsm) 


JANUARY 

STATION - 722696 


KM 

0.000 

1.000 

1.246 

2.000 

3.000 

4 .000 

5.000 

6.000 

7.000 
8 000 

9.000 
10.000 
11.000 
12.000 

13.000 

14.000 

15.000 


VAPOR P 
MEAN 
MB 

5.935 

4.616 

4.337 

3.462 

2.332 

1.351 

0.775 

0.461 

0.2S6 

0.130 

0.056 

0.026 

99.999 

99.999 

99.999 

99.999 

99.999 


SD VP 

MB 

3.818 
2.119 
1.830 
1.411 
1.162 
0.807 
0.511 
0.319 
0.175 
0.084 
0.039 
0.015 
99.999 
99.999 
99.999 
99.999 
99 999 


WHITE SANDS. NM 
SKEW VP TV 
MEAN 
DBG K 

1.82 281.09 

0.92 278.36 

0.70 277.47 

0.47 278.16 

0.73 273.00 

1.34 267.55 

1.81 261.25 

1.67 2S4.27 

1.31 246.79 

0.98 239.28 

0.78 231.77 

0.36 224.80 

999.99 219.39 

999.99 215.43 

999.99 213.97 

999.99 212.47 

999.99 210 13 


TV SKEW TV 
SD 

DEG K 

10.98 0.30 

7.64 0.14 

6.91 0.09 

5.31 -0.98 

5.06 -0.87 

4.76 0.42 

4.58 -0.32 

4.42 0.31 

4.36 -0.29 

4.23 -0.17 

3.80 0.01 

3.24 0.01 

3.16 0.32 

4.76 0.47 

4.99 -0.05 

4.03 -0.35 

3.58 0.33 


DEWPT T 
MEAN 
DEC K 
270.42 
268.06 
267 40 
264.56 
259.16 
252.32 
24 S . 97 
240.36 
234.44 
228.00 

219.99 
214.73 

999.99 
999.99 
999 99 
999.99 
999 99 


SD DPT 
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Perturbation Model 


GRAM data is based on monthly averages at a selected spatial scale. 
Variability at smaller scales is accomplished by the perturbation 
model 

The observed variability is partitioned into a large-scale (e g. large 
weather systems with time scales of several days) and a small-scale 
(e g. storms and turbulence) 

The large-scale is simulated with a cosine model to represent the 
wave nature of this phenomenon using a randomized phase 

The small-scale is simulated as a stochastic (random) process using 
a one step Markov technique 

- NewValue = OldValue * Correlation ♦ RandomFunction 

- Correlation decays exponentially with time and distance 

Monte Carlo runs of GRAM reproduce the observed monthly means 
and standard deviations 
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GRAM Produces Dryden Power Spectrum 



“Patchy” Parameter 


• Enables severe turbulence for patchy = 1 and only light- 
to-moderate turbulence otherwise (e g. patchy = 0) 

• For severe turbulence, the variances are increased by a 
factor of 6 to 12 (sigmas are increased by factor or 2.5 to 
3.5) depending on height 

• The probability of encountering severe turbulence in 
GRAM is consistent with the likelihood observed in 
nature (~ 0.2 to 2.5%, depending on height) 
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Patchy Turbulence 



Getting Started 


Contents of the GRAM99 CD 

• Reauests for GRAM99 


• OrderGUACA.txt 

- Sopo Yung at 256-544-0900 


- How to order full GUACA data CD 

- sopo.yung@msfc.nasa.gov 


• GRAMfixVer3.doc 

- Description of updates from earlier GRAM-99 

• Written aareements 


• README.txt 

- NASA employees: Software Request Form 


- Overview of CD contents 

- Domestic non-NASA employees: Software Usage 


• PCFiles directory 

Agreement (SUA) 


- Source code, PC executables, data files, README and other 
documentation, reference input/output test files 

- Foreign requests: SUA Form F (Export Controlled) 


• UnxFiles directory 

- Forms available from GRAM Web Site (see Chart 6) 


- UNIX-speafic source code and reference input/output test files 

7«H007 S3 


7/300007 


4 


PCFiles Directory 


PCFiles Directory (cont’d) 

• NameVer3.txt - Example input file 


■ READMEO.txt - General information as well as a description of the 

• OutVer3.txt - Resulting standard output file 


following readme files 

• SpclVer3.txt - Resulting special output file 


• README1 txt - Instructions for compiling and running the program 

• SpcsVer3.txt - Resulting species output file 


• README2.txt - Description of the NAMELIST data input file 

• guaca txt - Description of the Global Upper Air Climatic 


• README3.txt - Parameters available for the "Special Output" file 

Atlas data 


• README4 txt - New Global Gridded Upper Air Statistics (0-27km) 

• gguas.txt - Description of Global Gridded Upper Air 


data option 

Statistics files 


• README5 txt - New Range Reference Atmosphere Data option 

• gramfix.txt - Description of any fixes made since original 
code release 


and data files 

• README6 txt - New MET-99 model and wave perturbation model 
features 

• gramhist.txt - Development history 

• README.txt files 


• README7.txt - New option for user-selected initial perturbations 

7/202007 25 ^ 
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PCFiles Directory (cont’d) 

• GRAM95.pdf - GRAM-95 Tech Memo 

• GRAM99.pdf - GRAM-99 Tech Memo 

• GRAM99.exe - GRAM-99 PC executable 

• GRAMTRAJ.exe - Example PC executable for GRAM- 
99 incorporated into trajectory program 

• bldtraj.exe - Program to build array of GRAM input 
times/locations, as “trajectory" file input 

• atmosdat.txt - Atmospheric input file for GRAM 

• randin.txt - Sample random seed input file 

• Range Reference Atmosphere Data (1983) 

• FORTRAN source code ( f extension) 
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PCFiles Directory - GRAM99 Source Code Files 

D 

gram99.f 

- GRAM-99 mam program 

2) 

gramsubs.f 

- Various subroutines used by GRAM 

3) 

guaca. f 

- Subroutines related to GUACA data 

4) 

initial.f 

- Subroutine to initialize data 

5) 

met99prg.f 

- Source code for MET model 

6) 

models.f 

• Various sub-models use in GRAM 

7) 

random.f 

- Subroutines for random number generator 

8) 

rramods.f 

- RRA data handling routines for GRAM 

9) 

speconc.f 

- Subroutines for species concentrations 

10 ) 

bldtraj.f 

- Auxiliary program to build * * trajectory ” input file 

11 ) 

gguasrd.f 

- Auxiliary program to read and convert GGUAS data 

12 ) 

gramtraj.f 

- Routine to facilitate calling GRAM from within user’s 


driver program 
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PCFiles Directory - RRA Data Files 

• rrasites.txt - List of RRA sites, and 3-character 
site codes 

• tablel .xxx - Wind data files for RRA site xxx 

• table2.xxx - Pressure, density, temperature data 
files for RFtA site xxx 

• table3.xxx - Moisture data files for RRA site xxx 


PCFiles Directory - GUACA Binary Data Files 


GUACA monthly mean data files mxxxOI dat 
GUACA monthly standard deviation files sxxxOI dat 
Separate GUACA Files for Parameters xxx = 

- den = Density 

- dwp = Dewpoint temperature 

- hgt = Geopotential heights for pressure levels 

- sip = Sea-level pressure 

- tmp = Temperature 

- uwd = Eastward wind component 

- vwd = Northward wind component 

Data provided for January (Month 01) Penod-of-Record only. See 
file OrderGUACA txt for now to order full GUACA data CD 

See file guaca txt for how to set up GUACA data files in appropnate 
directory structure 
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GRAM Output Files 


Species Output File: 

Height 


Ocm> 

Lati- 

Long. 

con an 

Numr-er 


Concen- 

Number 



Tiim 

tude 

(E*w -1 

tta'- ion 

Density 


tration 

Density 

species 


(deg) 

(deg) 

-ppoiv. 

'* ra ** i ’ 


(ppmv) 

(«/m**3) 



140.00 

.470 

-164 .550 

0 OoOB*ilO 

0.0008*00 


0.0008*00 

0.0008+00 

H20 

1 




y. uf»uE*oo 

0 0008*00 


0.0008*00 

0.0008+00 

N20 

i CO 




. no 3E»00 

0.0008*00 


0.0008*00 

0.0008+00 

CH4 

j C02 




c Si ’8*05 

r ■ . 7278*1* 


9.1888*04 

9 . 4 70E+15 

N2 

i 02 




^ l . 

-• uO&B+i* 


2.3148*03 

2.38S8+14 


1 Ar 




1 

l . 9648*1 > 


9.7028-06 

1.0008+06 


i » 

139.80 

.510 

-164.430 

fl nunE*00 

0.0008*00 


0.0008*00 

0 . OOOE+OO 

H20 

1 03 




0 0iT'E+00 

0 . 0008*00 


0.0008*00 

0.0008*00 

N3D 

i co 

6 . 0 



" 00'iE*00 

0. 0008-00 


0.0008*00 

0.0008+00 

CH4 

j CO 2 




o 53 8*05 

* . ei7E*io 


9. 208E+04 

9.6098+15 

N2 

i 03 




. .5 *.-8*05 

t 5 >28*16 


2 . 3258+03 

2.4268+14 

0 

j Ar 




1 8908»02 

l .9728*13 


9. 5838-06 

1.0008+06 

He 

i H 



Output Files (cont’d) 


Standard Output File: 

Height Lati- Long 

(km) tude IB+W-) 

fine (deg) (deg) 

icl /RPA /Weight 

Pressure 

/Vap Pi 

(Ht/m—2) 

Tempera - 
Density/ ture/ 

vap Dene, Pewpt . 

(kg/m**3) (K) 

K-W 

wind 

(m/s) 

N-S 

Wind 

(m/s) 

Vert. 

wind 

(m/a) 


140.00 

.470 

-164.550 

9.8208-04 

4.3418-09 

690.1 

7.7 

.3 

.000 

Mean 




36. 3* 

13.3% 

23.3% 




M- 76 

.r+ 

OKK 

.0600 

.1% 

3.9% 

-6.9% 

83.7 

-9.4 


ranS 




2.9% 

6.2% 

12.6% 

45.0 

45.0 


sigS 




6.1% 

4.9% 

.7% 

-31.2 

-24.2 


ranL 




12.2% 

3.7% 

7.4% 

52.3 

52.3 


sigL 




6.2% 

8.8% 

-6.2% 

52.5 

-33.7 

7.27 

ranT 




12.5% 

7.2% 

14 .6% 

69.0 

69.0 

8.33 

s igT 




1.0438-03 

4.7228-09 

647.2 

60.2 

-33.4 

7.27 

Tot. 




44 . 8% 

23.2% 

15.6% 




T-76 




o.oonK+jn 

0 . 0O0E+00 

.0 



.0% 

H.Vj 




0.0008*00 

0.0008*00 

.0 



09 

pkih 



Output Files (cont’d) 

Special Output File: 


Tlae Hgtkm 
.00 140 000 

6.00 139.800 

12.00 139 600 

18.00 139 400 

24.00 139 200 

30.00 139.000 

36.00 138 800 

42.00 138 600 

48.00 138.400 

54.00 138.200 

60.00 138.000 

66.00 137 800 

72.00 137 600 

78.00 137 400 

84.00 137 200 

90.00 137 000 


Latitude Longitude DensMean PresMean 


47000 -164 
.51000 -164 
.S5000 -164 
.59000 -164 
.63000 -164 
.67000 -163 
.71000 -163 
.75000 -163 
.79000 -163 
.83000 -163 
.87000 -163 
.91000 -163 
.95000 -163 
.99000 -162 
1.03000 -162 
1.07000 -162 


.55000 4 .3410B- 
.43000 4.39698- 
.31000 4 . 45398- 
.19000 4.S119B- 
.07000 4.S710B- 
.95000 4.63118- 
.83000 4.69238 
.71000 4.75478- 
.59000 4.81828- 
.47000 4.88298- 
.35010 4.94918- 
.23010 5 . 01628- 
.11010 5 . 08478- 
.99010 5. 15448- 
.87010 5.22558- 
.75010 S . 29798- 


09 9.81988-04 
09 9.90218-04 
09 9.98558-04 
09 1.00708-03 
09 1.01568-03 
09 1.02428-03 
09 1.03308-03 
09 1.04198-03 
09 1.05098-03 
09 1.06018-03 
09 1.06948-03 
09 1.07888-03 
09 1.08838-03 
09 1.09798-03 
09 1.10778-03 
09 1.11768-03 


Tmean 

690.07 

687.28 
684 .48 
681.68 
678.87 
676.05 
673.22 
670.39 
667.55 
664.71 
661.89 
659.03 
656.16 

653.29 
650.42 
647.S4 


EWmean 

7.71 
7.50 
7.28 
7. OS 
6.84 
6.62 
6.39 
6.15 

5.90 
5.65 
5.42 
5.17 

4.90 
4 .63 
4.36 
4 .09 
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Controlling Special Output 


• Specify variables for output in models, f 

Write Uopp, 9000) elt.h, phi. thet, dgh, pgh, tgh.ugh, vgh.dh.ph, th.uh. ATMDSSSc 

k vh.drh.prh. sdh. sph. 0 . 01*ath*tgh, suh, svh.wh, leev ATMDS56 

9000 Format (P10.2,P9.3,P10.5 ) Pll.S,lp,2811.4,Op < P7.2,2P8.2,lp, ATMDSS7 

k 2811 . 4 , Op, F7 . 2 , 2P8 . 2 , 8P7 .2,13) ATMDSSB 

• Header can be changed in initial.f 

Write Uopp. 954) INIT208 

954 Formate Time Hgtkm Latitude Longitude DensMean', INIT209 

k PresMean Tmean EWmean NSmean DensPert PresPert INIT210 

k • Tpert BWpert NSpert Dpert% Ppert% SDdent SDprst • , INIT210a 

k ' SDtemK SDuwnd SDvwnd Wpert SeV) INIT211 



Sample Variables For Special Output 


Input Files 

• elt- time, sec 


• NAMELIST input file (mandatory) 

• h - height, km 


• GUACA database (full set must be ordered) 

• ph - total pressure, N/m**2 



• dh - total density, kg/m**3 


• Atmosdat.txt file (provided) 

• th - total temperature, K 


• Ranqe Reference Atmosphere input file (Data 

• uh - total east wind, m/s 


supplied, use optional) 

• vh - total north wind, m/s 


• Traiectorv input file (optional) 

• wh - total vertical wind, m/s 

• Complete list in file README3.txt 


• Random seed input file (Monte Carlo option) 
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Earth-GRAM99 Input File 


$ name in. da t 

atmpath - 'D:\GRAM99\atmosdat.txt' 

! path name Cor atmospheric data file 

guapath » 'O:\CDguaca' 

! path name for GUACA or GGUAS files 

trapath - 'null' 

! path name for trajectory input file 

prtpath - 'OutVer3.txt' 

('null' if none) 

t path name for standard formatted 


nprpath - 'SpclVer3.txt' 

output file ('null' if none) 

! path name for the ■special'' format 

conpa t h « ' SpcsVe r3 . txt • 

output file ('null' if none) 

! path name for species concentration 

rndpath • 'null* 

output file (‘null 1 if none) 

! path name for file containing more 

rrapath - 'D:\GRAH99\RRAdata' 

random number seeds (if needed) 
! DIRECTORY for RRA data 


hi - 140. 

t initial height (km) 


phil - 0.47 

! initial latitude (deg, North ♦ ) 


checi - -164.55 

1 initial longitude (deg, East ♦ ) 


flO - 230. 

1 10.7-cm flux 


flOb - 230. 

1 mean 10.7-cm flux 
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Earth-GRAM99 Inputs (cont’d) 


ap > 20.3 

t geomagnetic index 


mn » 1 

! month (1-12) 


ida - l 

t day of month 


iyr - 2006 

! 4-digit year, or 2-digit year yy, with >56>19yy 
«57-20yy 


ihro - 0 

l initial UTC (Greenwich) time hour (0-23) 


mino » 0 

l initial UTC (Greenwich) time minutes (0-59) 


seco - 0.0 

! initial UTC (Greenwich) time seconds (0.0-60.0) 


dphi - 0.04 

1 latitude increment (deg, Northward positive) 


dthet >0.12 ! longitude increment (deg. Eastward positive) 


dhgt - -0.2 

! height increment (km, upward positive) 


nmax > 701 

t maximum number of positions (including initial 



position,- ignored for trajectory input (iopt > 

0)) 

delt - 6.0 

! time increment between positions (real seconds) 


iopt - 0 

! trajectory option (0-no trajectory data,- otherwise 
unit number for trajectory input file) 

iopp - 17 

! "special* output option (0-no "special" output; 



otherwise unit number of "special" output file) 

iuO - 0 

! unit number for screen output (normally 6 or 0) 
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Earth-GRAM99 Inputs (cont’d) 

iup - 6 

! unit number for standard formatted output file 
(0 for none) 

ius - 3 

! unit number for atmosdat data 

iuc - 4 

! unit number for concentrations output (0 for none) 

iug - 22 

! unit for GUACA or GGUAS input data, 0-27km (0 for no 
GUACA or GGUAS data) 

iguayr - 1 

! Use: 1 for GUACA period of record, 2 for actual GUACA 

year (1985-1991), based on input value of year (iyr) 
3 for binary data converted from ASCII Global 
Gridded Upper Air statistics (GGUAS) POR data 
(conversion done with GGUASRD program, provided) 

iopr - 1 

! random output option (1 - random output, 2 - none) 

nrl - 1234 

1 first starting random number (1 to 9 • 10" *8) 

iun - 0 

l unit number for more starting random numbers (0 for none) 

rpscale - 1, 

.0 ! random perturbation scale, nominal-1.0, max-2.0, min-0.1 

iurra - 0 

! unit number for Range Reference Atmosphere (RRA) data 
(0 if none used) 
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Earth-GRAM99 Inputs (cont’d) 

sitelim • 

> 2 . S 

! lat-lon radius (deg) from RRA site, outside which RRA 
data are NOT used 

sitenear 

- 0.5 

l lat-lon radius (deg) from RRA site, inside which RRA 
data is used with full weight of 1 (smooth transition 
of weight factor from l to 0 between sitenear and 
sitelim) 

initpert 

- 1 

t Use 1 for user-selected initial perturbations or 0 
(default) for GRAM-derived, random initial 
perturbation values 

rpinit - 

10.0 

1 initial pressure perturbation value (% of mean) 

rdinit > 

1S.0 

! initial density perturbation value (% of mean) 

rtinit - 

-5.0 

! initial temperature perturbation value (* of mean) 

ruinit - 

3. 

! initial eastward velocity perturbation (m/s) 

rvinit - 

5. 

t initial northward velocity perturbation (m/s) 

rwinit - 

1 . 

i initial upward velocity perturbation (m/s) 

patchy - 
$End 

1 . 

! not equal 0 for patchiness,- 0 to suppress 
patchiness in perturbation model 
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Laptop Demo 2 - Horizontal Flyover 


Trajectory Input File 

Placeholder 


• Used for a known trajectory or an array of input positions 

• Overrides times and positions in the NAMELIST input file. 

• Set nmax equal to 0 in the NAMELIST input file. 

• Set iopt to a unit number in the NAMELIST input file. 

• Set trapath to the path name for the trajectory input file. 

• The trajectory input file consists of a text file (with no header 
record) containing columns of 

- time (sec) 

- height (km) 

- lat (+/- 90°) 

- Ion (+/-360 0 ), e.g. 
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Trajectory Input File Example 


0 . 00000E+00 

2.7854 IE- 02 

2 . 86084E+01 

-8 . 06041E+01 

2 . 02000E+ 00 

3 . 34559E-02 

2 . 86084E+01 

-8 . 06041E+01 

3 . 46000E+00 

5.23 94 7E- 02 

2.86084E+01 

-8 . 0604 1E+01 

4 . 54 000E + 00 

7 . 51876E-02 

2 . 86084E+01 

- 8 . 06041E+01 

5 . 44000E+00 

9 . 99682E- 02 

2 . 86084E+01 

-8 . 0604 1E+01 

6 . 22000E+00 

1 . 25826E- 01 

2 . 86084E+01 

- 8 . 06 04 1E+01 

6 . 90000E+00 

1 .51748E-01 

2 . 86084E+01 

-8 . 06041E+01 

7 . 52000E+00 

1 . 7814 8E- 01 

2 . 86 084E+ 01 

-8 . 06041E+01 

8 . 10000E+00 

2 . 05251E-01 

2 . 86084E+01 

-8 . 06041E+01 

8 . 64 000E+ 00 

2 . 32591E-01 

2 . 86084E+01 

-8 . 06041E+01 

9 . 16000E+00 

2.60852E-01 

2 . 86084E+01 

-8.0604 1E+01 

9 . 64 000E+ 00 

2 . 88638E-01 

2 . 86084E+01 

-8 . 06041E+01 

1 . 01000E+01 

3 .16808E-01 

2 . 86084E+01 

-8 . 06040E+01 

1 . 05400E+01 

3 . 4 5175E- 01 

2 . 86085E+01 

-8 . 06040E+01 
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Laptop Demo 3 - trajectory input 


• Placeholder 



Random Seed Input File 

• Used for Monte Carlo runs 

• A text file with a single column of integers that need not 
be random but should not repeat numbers 

• Set iun to a unit number in the NAMELIST input file. 

• Set iopr = 1 on the NAMELIST input file 

• Set rndpath to the path of the random seed file in the 
NAMELIST input file. 
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Random Seed Input File Example 


• The first number to be used is the value of nrl 
on the NAMELIST input file 

• The following file would produce 10 Monte Carlo 
runs: 

2 

3 

4 

5 

6 

7 

8 

9 

10 
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Laptop Demo 4 - Monte Carlo run 


Placeholder 


Compiling GRAM99 as a ‘Stand-Alone’ Program 

• On a PC. if you don't change anything in the source code, you don't 
have to re-compile, you can use the PC executable (gram99 exe) 
provided on the GRAM CD 

• Example of command-line compiling onaP£ with compile-and-link 
command “fl32“ 

- fl32 gram99 f gramsubs f guaca f initial f met99prg f models f 
random f rramods f speconcf 

- Creates PC-executable gram99 exe 

• Example of command-line compiling under UNIX with FORTRAN 77 
command *f77* 

- f77 aram99.f gramsubs.f guaca f initial f met99prg f models.f 
ranoom.f rramods. f speconc f 

- mv a. out gram99.x 

- Creates UNIX-executable gram99.x 


Running Stand-Alone GRAM 


Incorporating GRAM99 Into A Trajectory Code 

• On a PC (using executable gram99.exe) 

- Open a command-line window (Start/Run cmd.exe) 

- Type in cd GRAMdirectory to move to directory where 
executable resides and then type in gram99.exe, OR 

- Type in pathname\gram99.exe, where pathname is directory 
where executable resides 

- Type in the name of the desired NAMELIST input file (must be in 
same directory) 

• Under UNIX (using executable gram99.x 

- Type in cd GRAMdirectory to move to directory where 
executable resides and then type in gram99 jt, OR 

- Type in pathname\gram99.x, where pathname is directory 
where executable resides 

- Type in the name of the desired NAMELIST input file (must be in 
same directory) 


• Subroutine gramtraj provides a GRAM interface, 
accessible from within a user’s program by a “Call" 
statement 

• Subroutine gramtraj, along with other GRAM subroutines 
(gramsubs.f, guaca. f, initial. f. met99prq.f, models.f, 
random. f, rramods.f, speconc.f) should be compiled 
and linked with the user’s program 

• File gramtraj. f contains the gramtraj subroutine code, as 
well as some “dummy" code illustrating how a user 
program can be set up to call GRAM 

• File gramtraj. f has many comment statements, designed 
to aid the user in setting up the gramtraj subroutine in his 
or her program 
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Operating Tips - Verify Operation 
with Sample Input/Output Provided 


Operating Tips - Error Messages 

• After GRAM is properly configured, run it 
(as in chart 51) with example NAMELIST 
input file NameVer3.txt 

• If you get errors, see additional tips in 
following charts 

• If GRAM runs error free, verify output by 
comparing with sample output file 
NameVer3.txt, provided on GRAM CD 


• GRAM writes out error messages when it 
encounters certain kinds of error conditions 

- Use a text editor to try to find where in the code this 
error message is coming from; this may help in 
figuring out and correcting what caused the error 
condition 

- Extensive COMMENT statements are embedded 
throughout the code; these also may help in figuring 
out what the program is doing where the error 
occurred, and how to correct the error condition 

- If you wish to ask MSFC for assistance in resolving an 
error, it is VERY helpful to provide the input file you 
used, and the output you got (or the specific error 
message that you got) 
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Operating Tips - GRAM Perturbations 


Operating Tips - GUACA Data Problems 


• GRAM uses a one-step process to compute perturbations 

- If you change the random seed value, you will get a different sequence of 
perturbations 


• Read file guaca.txt (on the GRAM CD) 

• Verify that you have put GUACA files in the right pathname structure; 
examples: 


- If you repeat a given random seed value AND use the same set of calling step 
sizes, you will get the same sequence of perturbations 

- If you repeat a given random seed value, BUT you change the set of calling 


- GUACA data files in directory D:\GUACA\por\01 for January POR data 
on a PC and input guapath = D:\GUACA' 

- GUACA data files in directory /username/GUACA/por/01 for January 
POR data on a UNIX platform and input guapath = '/usemame/GUACA' 

• Possible “big-endian vs. little-endian" or other problems with binary data 

- Several GUACA read parameters can be changed in the code, and 
MAY solve the problem (or some compilers may have switches that 
correct a “big-endian vs. little-endian" problem) 

step sizes, you will get a different sequence of perturbations (because the 
steps are different) 

• Theoretically, GRAM perturbations can be evaluated at any calling 
frequency (step size) 

- If large steps are used, successive perturbations will be essentially 
uncorrelated 


- The small-scale Dryden spectrum means that perturbation energy will be 
present at all frequencies, with diminishing energy as the steps get smaller 


- Parameters provided in code in directory PCFiles should work on a PC 
running Windows (with Linux, there may be a problem) 

(calling frequency increases) 

- It is counterproductive (and may cause spurious results, depending on how 
your program treats high-frequency dynamics) to call GRAM perturbations too 
frequently (it is recommended that successive step sizes be used that are no 
smaller than the size of the vehicle) 


- Parameters provided in code in directory UnxFiles should work on most 
UNIX platforms 

- See comments in README0.txt and in the code on following charts if 
you need to change GUACA read parameters 

™«°° 7 K 
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Code to Change GUACA Read Parameters 
(PC Version Shown) 


C. . . Sec Che following paramecer values according Co your syscea GRAM 

C characce rise ics GRAM 

lswap » o gram 

lblwd « 0 GRAM 

irlbw « 0 GRAM 

nhdr - 45 GRAM 

lur « 5 CRAM 

dirsep = *v GRAM 

endsep « 'V gram 

cerachar > • • GRAM 

aerscac = -screech 1 GRAM 

sysfora = -binary- GRAM 

C... open che NAMELIST inpuc file GRAM 

Write (•, •) ’ Encer NAMELIST input file name' GRAM 

Read ( * , 5 i naaef lie GRAM 

5 Format (A 1 GRAM 

Open(unic<iur.file°nanieflle,scacua>-old' . loseae-loerr) gram 

If (loerr .ne.OIThen GRAM 

Wrice ( *. *) " Error opening NAMELIST file 1 GRAM 

SCop GRAM 

Endlf GRAM 

C GRAM 


67 

(8 

69 

70 

71 

72 
71 

74 

75 

76 

77 

78 
78a 
78b 
78C 
78d 
78e 
7Bf 
78g 
78h 
781 

79 
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GUACA Read Parameters (cont’d) 


C. . . Common variable lur is Che unit number of Che NAMELIST lnpue file GRAM 80 
C.... Common variable lswap is a swicch to control byte-swapping in the GRAM 82 
C GUACA data, which nay be required for certain platforms, e.g. GRAM 83 
C swapping is required for SGI machines. Syscea level commands may GRAM 84 
C also be used Co do Che byte swapping on some machines. To prevent GRAM 85 
C swapping, (e.g. on a VAX) set iswap»0. GRAM 86 
C.. Common variable lblwd is the block length, in 4 -byte words, if GRAM 87 
C the GUACA files have to be pre-read into a scratch file (lblwd * GRAM es 
C 0 means no pre-reading necessary; on a VAX use lblwd ■ 128) GRAM 89 
C.... Common variable irlbw is 4 If unformatted file record lengths are GRAM 90 
C specified in bytes (e.g Sun, RS/6000) ; irlbw la 1 if unformatted GRAM 91 
C file record lengths are specified in 4-byte words (e.g. SGI. VAX),- GRAM 92 
C Compiler options (e.g. old_rl on the SGI) can also change this. GRAM 93 
C... If GUACA files are to be read as direct access, fixed record GRAM 94 
C length, use lblwd record length in 4 -byte words and irlbw * l or GRAM 95 
C 4. If files are to be read as sequential access, use irlbw « o GRAM 96 
C and iblwd = 0. GRAM 97 
C. . . . Common variables dirsep, endsep and termchar are Interior and GRAM 98 
C ending directory separator characters and terminating character GRAM 99 
C in file pathname structures (e.g. use the following values for GRAM100 
C the various systems and file pathname structures) GRAM101 
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Earth-GRAM07 


• Currently undergoing Beta-testing 

• Expected release date November 2007 

• Operationally similar to GRAM99 
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Earth-GRAM07 New Features 


1. Revised Range Reference Atmosphere 
(RRA) data 

- In 2006, the Air Force Combat Climatology Center 
(AFCCC) developed a set of revised Range 
Reference Atmosphere (RRA) data including several 
new sites. 

- Earth GRAM-07 has the option of using either the 
2006 revised RRA data, or the earlier (1983) RRA 
data, as a replacement for conventional Earth 
GRAM climatology. 



GRAM07 RRA Sites (1983) 


• Ascension Island. Atlantic 

• Barking Sands. Hawaii 

• Cape Canaveral, Florida 

• Dugway Proving Ground (Salt Lake City). UT 

• Edwards Air Force Base, California 

• Eglin AFB. Florida 

• Kwajalein Missile Range, Pacific 

• Point Mugu Naval Air Weapons Center, CA 
■ Taguac. Guam 

• Vandenberg AFB. California 

• Wallops Island. Virginia 

• White Sands. New Mexico 

• Fairbanks. Alaska 

• Nellis AFB. Nevada 

• Shemya, Alaska 

• Thule, Greenland 

• Wake Island, Pacifio 

• Kodiak. AK (unofficial Developed by MSFC) 



GRAM07 RRA Sites (2006) 


Argentia Newfoundland (Sf Johns Airport) 

• Ascension Island. Atlantic 

• Barking Sands. Hawaii (Lihue) 

• Cape Canaveral. Florida 

• China Lake Naval Air Weapons Center CA 

• Dugway Proving Ground (Salt Lake City). UT 

• Edwards Air Force Base. California 

• Eglin AFB. Florida 

• El Paso. Texas 
Fairbanks. Alaska 

Huachuca Elec Prvng Gmd (Tucson). A Z 

• Great Falls MT 

Kwajalem Missile Range. Pacific 

• Nimes-Courbessac. France (STS TAL Site) 
Nellis AFB. Nevada (Mercury) 

• Point Mugu Naval Ax Weapons Center. CA 

• Roosevelt Roads (San Juan). Puerto Rico 

• Taguac. Guam (Anderson AFB) 

• Vandenberg AFB. California 

• Wallops Island. Vxgxtia (NASA) 

• White Sands Missile Range. New Mexico 
Yuma Proving Ground. AZ (San Diego. CA) 



New Features (cont’d) 


New Features (cont’d) 

2. Optional auxiliary profile input 


3. Updated thermosphere models 

- In addition to RRA options, an "auxiliary profile" 


- Users now have the choice of 3 thermosphere 

feature has been implemented. 


models: 

- This allows the user to input a data profile of 


- The revised Marshall Engineering Thermosphere 

pressure, density, temperature, and/or winds versus 


(MET-2007) model. 

altitude, in place of conventional climatology values. 


- The Naval Research Labs Mass Spectrometer, 

- Mean conditions are given by the profile if the 


Incoherent Scatter Radar Extended Model for the 

desired point is within a prescribed radius of 


thermosphere (NRL MSIS E-00) and the associated 

influence and are otherwise given by Earth GRAM 


Harmonic Wind Model (HWM-93). 

climatology. 


- The Jacchia-Bowman 2006 thermosphere model 



(JB2006). 
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New Features (cont’d) 


New Features (cont’d) 

4. Perturbation model revisions 


5. Coordinate system changes and revised earth 

- A new feature to update atmospheric mean values without 


reference ellipsoid 

updating perturbation values. 


- Equatorial and polar Earth radii for the "sea-level" reference 

- The ability to simulate large-scale, partially-correlated 


ellipsoid have been updated to World Geodetic System (WGS 

perturbations as they progress over time for a few hours to a 


84) values used by the GPS navigation system 

few days. 


- These are also equivalent (to 10 significant figures) to the 

- A multiple-trajectory driver routine that allows multiple 


Geodetic Reference System (GRS 80) values 

trajectories and perturbations to be simulated in one run. 


- Input values of altitude greater than 6000 Km are treated as 

- A multiple-profile driver routine that allows multiple profiles and 


geocentric radius values, rather than heights 

perturbations to be simulated in one run, with small-scale 


- Although all input latitudes are geocentric. Earth GRAM-07 now 

correlations maintained between the profiles. 


gives both geocentric and geodetic values on the output file. 

- Modifications which now produce wind shears from the model 


- A new subroutine has also been added which computes 

more similar to those observed at KSC 


horizontal distance from great-circle distance between two 
input latitude-longitude positions. 
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Abbreviations 


BACKUP CHARTS 
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• AFB - Air Force Base 

• DSNE - Design Specification for Natural Environments 

• ECMWF - European Center for Medium-range Weather Forecasting 

• EAFB - Edwards Air Force Base 

• GPS - Global Positioning System 

• GRAM - Global Reference Atmospheric Model 

• GGUAS - Global Gridded Upper Air Statistics 

• GUACA - Global Upper Air Climatic Atlas 

• MAP - Middle Atmosphere Program 

• MET - Marshall Engineering Thermosphere 

• POR - Period-of-Record 

• RRA - Range Reference Atmosphere 

• UTTR - Utah Test and Training Range 
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Data Set Period of Record 


• GUACA data: 1980-1991 

• GGUAS data (option): 1980-1995 

• 1983 RRAs:~ 1957-1979 

• 2006 RRAs: 1990-2002 
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Small-scale model 

Auto-correlated variable (density): 






Cross-correlated variables (pres-dens, temp-dens) 

Pi=r>Px+ r „Pi +r q Q 


7/20/2007 72 



Dryden Spectrum 


According to Lumley and Panofsky "The Structure of 
Atmospheric Turbulence", the energy spectrum is given by 

/r (* )= ik e '“ r( ' ) " V/ 

Where r is the correlation coefficient given by 
r(l) = e' ,,L for /£ 0 
r(l) = e" L for I <, 0 



L 

1 + k 2 L 2 
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Patchy Turbulence 



— Patchy=1 
Patchy=0 
• Se^re 
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Genesis Entry Profiles 




August GRAM Density vs NCEP Ciimatology (1979-2000) 

Lat 41 N Lon 113W 


0.95 1.00 1.05 1.10 

Density/1976 US Standard 




40.8 
40.7 
40.6 
£ 40.5 
40.4 
% 40.3 
| 40 2 
40.1 
40 0 
39 9 


Target Location; (246.55 deg. 40.3167 deg) 

Final Landing Location: (246.4776 deg. 40 3683 deg) 


O 


UTTR Boundaries 


Final 
Landi 
Location 


Landing j Landing Ellipse 


OD s06015f 
Landing Ellip 
w/ E-2 hr Winds's 


ODs06015f 4 

Landing Ellipse J 




245 8 246.0 246.2 246.4 246 6 246 8 247 0 247.2 247 4 
East Longitude, deg 




